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Abstract. We describe the evaluation of the anomalous dimensions of twist-2 deep in¬ 
elastic light cone operators to 0{1/Nf) as a check on future perturbative calculations. In 
particular we present recent results for the singlet gluonic operator dimension in polarised 
and unpolarised scattering and give three loop predictions for the 0{1/Nf) gluonic eigen- 
operator. The section of the 3-loop 0{1/Nf) DGLAP splitting function proportional to 
the adjoint quadratic Casimir is also calculated for the singlet gluonic operators. 
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1 Introduction. 


Our modern view of hadronic matter is based upon the rich quantum held theory of 
Quantum Chromodynamics (QCD) (for a review see [Q]) in which we picture hadrons as 
bound states of valence quarks in a sea of vacuum excitation quark/antiquark pairs and 
gluons. Because QCD has the remarkable property of asymptotic freedom we can use 
high momentum transfer reactions such as deep inelastic scattering {lepton + nucleon 
—> lepton + hadrons) to compare experiment with the predictions of perturbative QCD 
(pQCD) and so rehne our ideas on the structure of strongly interacting particles. 

The formal QCD approach to deep inelastic scattering (DIS) allows us to predict struc¬ 
ture function behaviour using moment sum rules. Key to this method is the formalism 
of the light cone expansion (LCE) in which the non-local time ordered product of elec¬ 
tromagnetic quark currents appearing in the general DIS cross-section is expanded in a 
series of local, spin-n operators together with c-number Wilson co-efficients. The LCE 
is easily seen to be dominated by operators of lowest twist (r = operator dimension — 
operator spin). For QCD the lowest twist operators available are twist-2. It can be shown 
by using the anomalous dimensions associated with the renormalisation of insertions of 
these operators in 2-point Green functions, that we can determine the evolution of the 
moment sum rules with the large momentum transfer scale. 

The unpolarised twist-2 operators are [§, 
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and the polarised twist-2 operators are |^, 
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Here, i/j is the quark held, is the gluon held strength tensor, is the QCD covariant 
derivative. A* are the generators of SU{Nf), S stands for symmetrisation over Lorentz 
indices and NS denotes the non-singlet operator under SU{Nf) as opposed to the other 
singlet operators. 

An alternative, less abstract means of providing theoretical predictions for DIS is 
given by the parton model. Here the structure functions may be expressed in terms of 
parton distribution functions (PDFs). These give the probability that the struck nucleon 
constituent (parton) carries a particular fraction of the total nucleon momentum. The 
evolution of the PDFs are then governed by the DGLAP equation [^]. This contains 
the DGLAP ‘splitting functions’ which express the probabilities for the struck parton 
undergoing certain collinear decays as the energy scale changes. The DGLAP splitting 
functions can be explicitly obtained from the previously mentioned operator anomalous 
dimensions through an inverse Mellin transform. Solutions of the full DGLAP equation 
provide us with scale dependent PDFs which may than be used as inputs for other hard 
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scattering processes where the formal QCD approach is not necessarily applicable as well 
as giving us a deeper insight into hadron structure in the asymptotic regime. 

With facilities such as HERA currently opening up new kinematical regimes for unpo¬ 
larised and (soon) polarised DIS, there presently exists a need to perform a NNLO twist-2 
analysis to ensure a high precision examination of pQCD. This is especially important 
where the new order’s effects are expected to be seen experimentally (as in low x, high 

scattering at HERA). This programme requires the 2-loop Wilson co-efficients and 
3-loop anomalous dimensions of all the DIS twist-2 light cone operators with full spin de¬ 
pendence. Although the complete 2-loop operator dimensions have been known for many 
years now i I, 0, |], the full 3-loop computations have only become viable in recent 
years. The current state of play is that the 3-loop unpolarised singlet and non-singlet 
operator dimensions have been calculated for particular spins of the operator (n = 2, 4, 6, 
8 and n = 2, 4, 6, 8, 10 respectively) [|, 0. In addition, the 2-loop hnite parts required 
for the full 3-loop results have recently been calculated for all twist-2 unpolarised and 
polarised operators BUI- As one might expect, these are particularly difficult calcu¬ 
lations involving thousands of Feynman diagrams. It is clear that a check on the hnal 
full 3-loop results would be useful. One way of approaching this would be to calculate 
the dimensions using an alternative expansion parameter to the QCD coupling and check 
that there is agreement where overlap exists. We can do this by use of the 1/A^ expansion 
and a critical point approach. 


2 The Method. 


By applying the critical large Nf method developed in the series of papers IjT^ it is possible 
to obtain expressions for the twist-2 operator dimensions to 0{1/Nf). The method entails 
the analysis of operator insertions in QCD Green functions at a d-dimensional non-trivial 
renormalisation group hxed point. This hxed point may be found as a stable zero of the 
d-dimensional four-loop QCD /d-function |T3 and is located at 
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where Oc is the strong coupling constant at criticality, C 2 {R) and C 2 (G) are the funda¬ 
mental and adioint quadratic Casimirs respectively, tr(T“T^) = 1(5“^ for T“ the generators 
of SU{N^) and d = 4 - 2e. 

Several good things come out of using such an approach; 


1. Propagators take on a simple, dressed, scaling form due to the scaling properties of 
Green functions at a renormalisation group hxed point. 
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For example, the quark and gluon propagators are respectively |TB 
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where b is the covariant gauge parameter, d = 2/i, A and B are momentum indepen¬ 
dent amplitudes and the exponents a and (3 are defined using simple dimensional 
analysis of the massless QCD action as 


a = - I + \r] , j3 = I - rj - X (2-3) 


Here, the critical exponent r] is the quark anomalous dimension associated with the 
quark wave function renormalisation at the critical coupling and similarly y is the 
anomalous dimension of the quark-gluon vertex. 


2. Using the critical renormalisation group it can be shown that the form of a particular 
0{1/Nf) operator anomalous dimension at Oc may be calculated by evaluating the 
residues of first order poles (we use an analytic regularisation of the gluon dimension) 
in appropriate 0{1/Nf) two-point Green functions at Oc. These Green functions 
consist of an insertion of the operator into a QGD two-point function. (We also 
need to take the field renormalisations into account). 


3. Since, in the language of statistical mechanics, we are working at a critical point, we 
know that we have an essentially massless theory which enables us to use integration 
tricks such as uniqueness and conformal transformations to evaluate diagrams. 


4. In d-dimensions and to 0{1/Nf) there exists a universality equivalence between 

H 


QGD and the non-abelian Thirring model (NATM) at Oc 


This may be seen 


through the reproduction of the QGD triple and quartic gluon vertices by fermion 
loop integration in three and four point NATM Green functions at Oc using the above 
propagators. The upshot of this is that we can use the simpler NATM interactions 
in our calculations and dispense with the tricky three and four-point gluon vertices. 


5. The fact that we work with a fixed d-dimensional spacetime with an analytic rather 
than a dimensional regularisation seems to alleviate some of the problems caused 
when treating 75 in arbitrary dimensions 0 - This obviously becomes important 
in the calculation of the dimensions of the polarised operators. 


3 Comparing singlet operator dimensions at 0{1/Nf) 
with perturbation theory 

A slight complication in this calculation arises due to mixing of the singlet operators in 
both polarised and unpolarised scattering. Since the singlet light cone operators share 
the same quantum numbers and have equal canonical dimension in strictly four dimen¬ 
sions, they mix under renormalisation. This means that we have to consider a matrix of 
renormalisation constants for these operators 

(3-1) 
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Here the indices i,j = q,g and we are refering to the operators ( | 1 . 1 |) and ( p.. 2 D . The 
anomalous dimensions, 'jijia), are defined by 

7 p(«) = ( M""] ] (3.2) 

where 7 *^( 0 ) = P{a){d/da) InZ^ and P{a) is the renormalization group function governing 
the running of the QCD coupling constant a. The entries in 'jij{a) depend on the colour 
group parameters, Nf and n and since it is the l/Nf corrections that we are interested in, 
we define the explicit form of the entries as 

Iqqifl) = o-ifl + {a2iNf + 022 ) 0 ^ + + 0-32^^ + 0 - 33 ) 0 -^ + 0{a‘^) 

jgqia) = bia+{b2iNf+ b22)a‘^ + {bsiN^+ b32Nf+ b33)a^ + 

Iqgio) = CiNfa + C2Nfa^ + {c3iNj + C32Nf + 033 ) 0 ^ + O(a^) 

lgg{(^) = (diiA/+ ^ 12)0 + (d 2 ii^ + ^ 22 ) 0 ^ + (<^31.^^ + <^32.^ + + 0(a^)(3.3) 

where Nf = T{R)Nf and the coefficients aij, bij, Cij and depend on n and the colour 
group Casimirs. 


We note that this matrix of anomalous dimensions has eigenvalues 


X±{a) — ^{'Jqq + 'ygg) ± ^ {'Jqq ^ gg) ^-^'^qg^gq 
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Expanding in powers of a and retaining the same orders in 1/Nf with the definitions 
we find. 
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Here we can see that A+(a), A_(a) are dominated by contributions from the gluonic and 
fermionic operators respectively. 


When in our approach we consider the singlet sector operators in d-dimensions at Oc, 
we find that they no longer mix. (It is easy to see that the difference is 0(e).) This means 
that by calculating in d-dimensions, we are accessing the above mixing matrix eigenvalues 
in perturbation theory. This is borne out by explicit calculations. By evaluating the 
graphs required for 'jgg(ac), Jqqiac) we actually obtain A+(ac) and A_(ac) respectively, 
with the universality equivalence between QCD and NATM at 0{1/Nf) accounting for 
the contributions from the off-diagonal dimensions. 


4 Results 

. The 

polarised non-singlet operator result was published in [l^ together with the expressions 
for A_(ac) to 0(1/A^) and A+(ac) to 0(1) for polarised and unpolarised scattering. 


The unpolarised non-singlet operator dimension at 0(1/A^) was published in |18 
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We now give new results 10.0 for A+(ac) to 0{1/Nf) for the unpolarised and po¬ 
larised gluonic operator together with the 3-loop n-dependent eigenvalue predictions for 
comparison with perturbation theory (|3.5|) . The relevant graphs are given in Fig. 1 with 
a different operator insertion Feynman rule [|^ ^ used for the separate unpolarised and 
polarised cases. 









I 

I 








Fig. 1. Leading order diagrams for A+(ac). 

For the unpolarised gluonic singlet operator we have, 

A+4(ac) = — [8/i^n^ -|- -|- 8p^ -|- 2ij?n^ + — 24/i^n 
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— 8/i^n^ -f 240/i^n^ -|- 244/i^n -|- 316/i^ -|- 16/i^n"^ -|- 32/i^n^ 
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implying, 

&31C1 _ 64(n^ + n + 2)^(S'i(n))^C2(-R) 
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For the polarised gluonic singlet operator we have, 
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Throughout these calculations we used reduce and FORM ||^ to handle tedious 
amounts of algebra. The quantity r/° is dehned by 

(2p-l)(/i-2)r(2p) 
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and we have set A+(ac) = •^+*(®c)/-A^*- The hnite sum S'i(n) is given by S'i(n) = 

E”.i 1/i- 

These expressions agree exactly with all known perturbative results. This can be seen 
by putting p = 2 — e in iW^ and (|4.4|) and expanding in powers of e/Nj. An interesting 
feature of the mixing matrix ®3) is that the 0(1/A^) contribution to 'jgq{a) depends only 
on the Casimir C 2 {R)- This is evident from the one and two loop results for all n and 
the three loop results for n < 8. Assuming this to be true for all n at three loops we can 
see from 


and (|4.4|) that we can deduce the exact form of the coefficients ^31 in the 
C 2 {G) sector. This means that by using an inverse Mellin transform we can calculate the 
0{1/Nf) part of the three loop DGLAP function Pgg for both unpolarised and polarised 
cases. For the unpolarised splitting function we obtain, 
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Similarly for the polarised splitting function, 
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These results complete the programme to calculate the 0(1/A^) corrections for the 
anomalous dimensions of the twist -2 light cone operators. At present it seems that a 
continuation of this programme to include contributions of 0 ( 1 /A^ 2 ) pg viable. We 
conclude by noting that the results may also be useful in estimating the full three loop 
corrections to the operator dimensions by using asymptotic Fade approximant techniques 
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